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Abstract. We describe a renormalization group transformation that is related to the 
breakup of golden invariant tori in Hamiltonian systems with two degrees of freedom. 
This transformation applies to a large class of Hamiltonians, is conceptually simple, and 
allows for accurate numerical computations. In a numerical implementation, we find a 
nontrivial fixed point and determine the corresponding critical index and scaling. Our 
computed values for various universal constants are in good agreement with existing data 
for area-preserving maps. We also discuss the flow associated with the nontrivial fixed 
point. 
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1. Introduction 

The renormalization group transformation considered in this paper acts on analytic Hamil- 
tonians for systems of two degrees of freedom. It is designed to apply to the stability 
problem for invariant tori associated with frequency vectors uj = (lo\, L02), for which UJ1/0J2 
is a reduced quadratic irrational. For simplicity, we will focus from the beginning on the 
case of the golden mean. Our formalism can also be extended to higher dimensions, as in 
[4], but the results may be qualitatively different; see also [15-17]. 

We start with some simple facts about the golden mean $ = | + \\f§. The numbers 
d and — l/i? are solutions of a quadratic equation with integer coefficients, namely the 
equation for the eigenvalues of the matrix T = [®\]- The corresponding eigenvectors 
are uj = (l/#, 1) and O = (1,— 1/#), respectively. The matrix T can also be used to 
generate the continued fraction approximants dk = bk/a-k of More specifically, we have 
(a k ,b k ) =T fc (l,l), for fc = 0,1,2,.... 

The Hamiltonians considered here are assumed to have a rapidly convergent Fourier- 
Taylor series of the form 

H(q,p)= H v , a F v , a (q,p), F v , a (q,p) = (uj-p)^(Q-p)^e i ^, (1.1) 

for all q in IR 2 , and for p in some open neighborhood of the origin in R 2 . The sum in 
this equation ranges over all multi-indices (v, a) in / = Z 2 x N 2 , with N being the set 
of all non-negative integers. We will refer to p and q as the action and (lifted) angle 
variable, respectively. Recall that the flow defined by H is (q,p) = (V2H, — Vi-ff), where 
V j denotes the gradient with respect to the j-th variable. 

In order to motivate our analysis, let us first consider the integrable Hamiltonian 
Ho(q,p) = uj ■ p + ^(O • p) 2 . When restricted to a fixed level set Ho(q,p) = E, this 
Hamiltonian has an invariant torus O^o whose frequency vector is proportional (in fact 
equal) to uj. In addition, for every sufficiently large integer k, we can find a periodic orbit 
Ok whose frequency vector is a constant multiple of (ak,bk). The same can be observed 
(and proved) for analytic Hamiltonians that are sufficiently close to Hq. In fact, the orbits 
Ok accumulate at the torus Ooo . To give a more precise description, assume for simplicity 
that each of the orbits Ok contains a unique point (qk,Pk) with qk = 0. This holds e.g. for 
many Hamiltonians that satisfy H(q,p) = H(—q,p). Then pk approaches a limit poo as k 
tends to infinity, and the distance dk between pk and p^ tends to zero with an asymptotic 
ratio dk+i/dk — > $~ 2 . This should not bee too astonishing, as |#fc — i?| tends to zero with 
the same asymptotic ratio. What is more astonishing is what seems to happen in many 
one-parameter families (3 1— > Hq + /3F, as the strength of the perturbation is increased: By 
looking for orbits with a given degree of stability (e.g. residue 1), among the orbits Ok 
for sufficiently large k, one finds a convergent sequence (3 m , Pm+i, • • • in some parameter 
interval J, with the following property. The Hamiltonian Hq + /3F, with (3 e J, has an 
orbit with the given stability only if (3 = (3k for some k > m, and in this case, the orbit 
is unique and equal to Ok ■ Furthermore, the ratios (fik-i — /3k) /(/3k — Pk+i) converge 
to a limit 5, as k tends to infinity; and the ratios dk+i/dk, for the Hamiltonian with the 
parameter value (3^ = limfc^oo (3k, converge to a limit A different from # -2 . The critical 
index 5 and the critical scaling A appear to be independent of the family considered, within 
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a large class of Hamiltonians. The observed values are 6 = 1.627 . . . and A = 0.326 — A 
more detailed description of these phenomena can be found in References [1,2,10-13]. 

Such scaling and universality properties suggest the use of renormalization group 
methods. The idea is to find a transformation 7Z, acting on a space of Hamiltonians, which 
contracts differences that are irrelevant for the phenomena being considered. According to 
the description above, 7Z should therefore have two fixed points: One that is a prototype 
of a Hamiltonian with trivial asymptotic scaling, attracting all other Hamiltonians of this 
type, and another one that plays the same role for Hamiltonians with critical scaling. The 
fixed point equation 71(H) = H should of course imply self-similarity for the set of 
orbits of H. A transformation 9ft which has the potential of satisfying all of the necessary 
requirements (modulo the existence of a trivial unstable direction) has been described in 
[4]. It is of the following form: 

&(H) = -H o U H o T M - e , T^q, p) = [Tq, fiT^p) . (1.2) 
A* 

Uh is a canonical change of coordinates that will be defined later; in particular, Uh is 
the identity map if H(q,p) does not depend on q. The constants e = e(H), r = t(H), 
and n = n(H) are determined by three normalization conditions for H = 9ft(i7). In this 
paper, we choose these conditions to be H 00 = 0, i?o,(i,o) = 1) an d #0,(0,2) = \- If they 
cannot be satisfied, then dt(H) remains undefined. It is easy to check that the integrable 
Hamiltonian H (q,p) = 10 ■ p + |(f2 ■ p) 2 is a fixed point for 9ft, with e(Ho) = 0, t(H$) = 
and fjL(H ) = tf~ 3 . 

The first step H 1— > H o Uh in the definition of dt(H) is similar in spirit to the 
averaging operation used in KAM theory [18,19]. Its purpose is to eliminate irrelevant 
(non-resonant) degrees of freedom from H. The next step H 1— > H0T1 represents a "shift" 
of frequencies, in the sense that if H has a periodic orbit for a frequency vector u, then 
H o Ti has a (trivially) related orbit for the frequency vector T~ l u. Notice that Ti is a 
canonical transformation. In order to re-normalize the resulting Hamiltonian, we perform 
a rescaling H fi~ 1 H(., /x.) of the action variable p, a rescaling H tH of the energy (or 
time), and subtract a trivial constant e. We note that the shift and scalings are directly 
related to the observation of self-similarity near invariant w-tori. They have also been 
used in other renormalization transformation [1-8]. 

The change of coordinates Uh will be chosen within a restricted class of canonical 
transformations. Each transformation U : (q,p) 1— > (q + Q,p + P) in this class is associated 
with a generating function 4>(q,p) that is 27r-periodic in the components of q, and U = 
is defined implicitly by the equation 

Q(q,P) = V 2 (f>(q,p + P{q,p)), 
P(q,p) = -W 1 <f>(q,p + P(q,p)). 

One of the consequences of this restriction is that Uh cannot include p-translations like 
V s : (q,p) 1— > (q,p — sQ). As a result, 3ft has a trivial unstable direction. This could be 
avoided easily by including suitable translations in the definition of our renormalization 
group transformation. More specifically, we can define 71(H) = (JR.(H o V s )) o V_ s , with 
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s = s(H) determined in such a way that H = TZ(H) satisfies #0,(0,1) = 0- Our numerical 
result for dt, described in the next section, can be translated into corresponding results for 
1Z. They indicate that 1Z has a trivial attractive fixed point (namely Hq), and a nontrivial 
fixed point ii* whose stable manifold is of codimension one. 

Concerning the notion of relevant or irrelevant degrees of freedom for analytic Hamil- 
tonians, we note that even a critical o;-torus is smooth in the direction of the flow. This 
suggests that the most relevant perturbations are the ones that distort the cu-torus, and 
the geometry of nearby orbits, mainly in the directions transversal to the flow. If we 
assume that (or work in coordinates where) the average velocity on the w-torus of H 
is proportional to u, then these perturbations are given by the functions F v ^ a whose fre- 
quency vectors v are almost perpendicular to u>. They are usually referred to as "resonant" 
modes. The idea is to eliminate all but these resonant modes. Another way of seeing why 
this should be done is the following: Assume that, as a function of the angle variable, H 
is analytic in a strip given by |Im ui ■ q\ < r and |Im Q ■ q\ < R. Then H o T± is analytic in 
a similar strip, but with r and R replaced by r$ and R/fl, respectively. Thus, in order to 
avoid that dt shrinks the domain of analyticity, we have to regularize H in the direction 
(Im q) oc ui. One way of doing this is by eliminating the non-resonant modes of H. To be 
more precise, we will call F v ^ a a resonant mode if (u, a) belongs to the set 

1+ = {(u,a) G Z 2 x N 2 : \u ■ u\ < <r\\v\\ + k|H|} , (1.4) 

or a non-resonant mode if (is, a) belongs to the complement I~ of this set in Z 2 xN 2 . Here, 
cr and k are positive constants to be chosen later, and |.| denotes the i 1 norm on M 2 , that 
is, 1 1 a; 1 1 = + ja^l- Correspondingly, 

I ± H= H »,«F v , a (1.5) 

(i/,a)e/± 

will be referred to as the resonant (+) and non-resonant (— ) parts of the function H. 
Notice that the non-resonant part of a function that only depends on the action variable 
is zero. 

Given a Hamiltonian if, we now choose Uh to be the canonical transformation 
whose generating function is defined implicitly by the equation 

I" (# o Efy) = , 1+0 = 0, (1.6) 

provided that a (unique) solution exists. This definition oiUu is similar in spirit to the one 
used in [4], but it is considerably simpler. In particular, if H is analytic and sufficiently 
close to Hq (in an appropriate sense), then equation (1.6) can be solved by a relatively 
simple Newton type iteration [9]. Concerning our definition of resonant modes, we note 
that there are other possible choices that should lead to the same results. But if I + is too 
small (e.g. k = 0), then the equation (1.6) becomes singular, and if I + is chosen too large, 
then 3? will not be analyticity improving. 

In the remaining part of this paper, we present some numerical evidence for the 
existence of a nontrivial fixed point of 3? and discuss some of the expected properties of 
this fixed point. 
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2. Results 

In order to simplify the problem of finding a nontrivial fixed point for dt, we will restrict 
our search to a suitable subclass of Hamiltonians that is left invariant by 3?. In particular, 
we will only consider Hamiltonians that are even in q, or equivalently, invariant under the 
reflection V : (q,p) i— > (—q,p). The equations (1.3) and (1.6) will be restricted accordingly, 
to generating functions cj) that are odd in q, and to canonical transformations Uh = 
that commute with V. We note that these parity properties of <p and follow from H 
being even in q, if H is close to H . As a further restriction, we assume that 

H(q,p) =uj -p + h(q,Cl-p) = u> ■ p + ^2h u , n cos(v ■ q)(fl ■ p) n , (2.1) 

with h v ^ n G R. For the generating function of the canonical transformation Uh = 
we make the corresponding ansatz <f>(q,p) = <p(q,0> • p), where <p(q,z) is a sum of terms 
<^„,n sm(v • q)z n . Given such a function <p, we can turn the second part of (1.3) into an 
equation for a function Z = Q ■ P, that only depends on q and z = O • p and is even in q: 

Z(q, z) = -Q- Vi^(<z, z + Z(q, z)) . (2.2) 

From a solution of this equation we obtain 

U<t>(q,p) = (q + d 2 <fi(q,z + Z(q,z))Q,p- Vi^(g, z + z))) , z = Q-p, (2.3) 

where c^V? denotes the partial derivative of <p with respect to the second argument. It is 
easy to check that if H is of the form (2.1), then the same is true (formally) for H o . 

Within this restricted class of Hamiltonians, the transformation dt involves the follow- 
ing steps: Given H, we first solve the equation I~(H o U$) = for a generating function 
(j) = of the type described above. After composing H with the canonical transforma- 
tion Uh = Uef, , the new Hamiltonian is again of the form (2.1), and made up entirely of 
resonant modes. Some of these resonant modes are turned into non-resonant ones in the 
next step: H \— > H o T\ . After this step, we perform the scaling H \— > r / u _1 if(., fi.) — e 
described after (1.2). The resulting function 31(H) belongs again to our restricted class 
of Hamiltonians. We note that t(H) = $ for all Hamiltonians H in this class, i.e., the 
renormalization of time is trivial. 

In our numerical implementation of these steps, we approximate functions like h(q, z) 
or (j)(q, z) by finite linear combinations of modes cos(z/ • q)z n or sin(z/ • q)z n , respectively, 
and represent them by arrays of Fourier-Taylor coefficients. The cutoff is taken to be 
of the form + n < N, and all intermediate results are truncated to the same fixed 
"degree" N. We have written subroutines for the addition, differentiation, multiplication, 
and composition of such functions; and similarly for the inverses of these operations, using 
Newton's method when necessary. We also use Newton's method to solve (1.6). The 
projections are fixed once and for all, by choosing a = 2/5 and k = a 2 . 

In order to find a good approximation for the expected nontrivial fixed point of 
3?, we started with some initial guess, and iterated (the numerical implementation of) 
a contraction that has the same fixed point as 3?. The contraction was chosen of the 
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form KoP, with V a projection onto a codimension two manifold passing through , 
defined by the equation p(dt(H)) = p(H), for a suitably M 2 -valued function p. For the 
degrees N = 7,9,11,13,15, we iterated the contraction, until reaching a Hamiltonian 
Hi N) (q,p) = u> ■ p + hi N) (q, O ■ p) that satisfies the fixed point equation for 3?numeric > U P to 
an error less than 3 x 10 -15 in the norm 

\\h\\ =sup|fc I/|B |e 2 IHI. (2.4) 

At degree N = 7, the initial guess was obtained by a bisection procedure, as described 
e.g. in [6-8]. Then Hi 7) was used as a starting point at degree 9, and similarly for higher 
degrees. As a result of this iteration procedure, we get approximate values for the Fourier- 
Taylor coefficients of if* , and for the scaling parameter (i(H*). In particular, the computed 
coefficients satisfy ||/ii JV) || = 0.92 . . . for all N, indicating that H* has a decent domain of 
analyticity. Our results for n(Hi N) ) are listed in the second column of Table 1. The A = 15 
value for the energy scaling factor differs by less than 10 -7 from the value 7.0208826 
obtained from a renormalization group analysis of area-preserving maps [10-13]. 

The fact that 3?oP was observed to be a contraction, indicates that in the space 
of Hamiltonians considered, 3ft has (at most) two unstable directions at . In order to 
determine these directions, we computed the image of a two dimensional subspace under 
iteration of a numerical approximation to _D3ft(.£/*). From the action of this map on the 
limiting "invariant" subspace, we computed approximate values for the two expanding 
eigenvalues 5i and 5 2 of D$l(H*); see columns 3 and 4 of Table 1. The value of 5 2 obtained 
in [13] is 1.6279500. As we mentioned in the introduction, one of the eigenvalues 5j 
is associated with an irrelevant degree of freedom corresponding to a translation in the 
action variable. Indeed, our numerical results are compatible with the expected relation 
Sin = —1/$, involving the scaling parameter [i for the action variables. We note that the 
limiting discrepancy of about 2 x 10~ 7 is due to the fact that we used a fixed difference 
quotient for our computation, instead of writing a program for the derivative of 3ft. 
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7 


0.230376748 


1.1 x 10" 5 


1.66289337 


9 


0.230454722 


6.8 x 10" 7 


1.62803658 


11 


0.230460077 


1.9 x 10" 7 


1.62795826 


13 


0.230460180 


2.1 x 10" 7 


1.62795515 


15 


0.230460194 


2.1 x 10" 7 


1.62795473 



Table 1. Scaling /u, and eigenvalues 6j of D5R(H«) 



We do not expect Ddt(H^) to have additional eigenvalues of modulus > 1, when 
considered on a space of general Hamiltonians of the form (1.1). Support for this comes 
from earlier numerical work [5], where a renormalization group transformation similar to 
3ft was implemented, without the restrictions considered in this section, except for the 
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condition that H (q, p) be an even function of q. At the approximate nontrivial fixed point 
of this transformation, no eigenvalues of modulus > 1 were found, besides 5\ and 62 ■ 

Next, we consider the scaling transformation Sh = Uh ^n(H) which maps orbits for 
^St(H) to orbits for H. We will assume that 3? has a fixed point H* of the form (2.1), close 
to our Hamiltonians Hi N \ Since Sh* commutes with the reflection V and maps the energy 
surface H*(q,p) = E to the surface of of constant energy (E + with |///#| < 1, we 

expect that Sh* has a fixed point of the form (0,Poo)- Our numerical results indicate that 
this is the case, and that p^ ~ (0.05616862,0.01352928). Since Uh* Ti is a canonical 
transformation and commutes with F, the derivative of Sh* at the point (0,Poo) has to be 
of the form 

A T B = fA, (2.5) 

where A and B are 2x2 matrices. By comparing derivatives of 3? (if*) and at the 
point (0,poo), one also finds that A has two eigenvalues Ai = d and A2 , with eigenvectors 
v\ = V2-ff*(0, Poo) and v 2 = fi, respectively. Thus, since A T B = /xl, the eigenvalues of 
-B are A3 = and A4 = fi/fl, and the corresponding eigenvectors are V3 and V4 = u>, 

with w 3 perpendicular to v± . We note that the triviality of Ai and V4 is related to the 
fact that the geometry in the direction of the flow is trivial, and that the vector field for 
H* is invariant under translations (q,p) 1— > (q,p — soj). The approximate values for A2 , 
obtained from our approximate fixed points H* , are listed in column 3 of Table 2. At 
= 15, our value for I/A2 , when truncated to 8 digits, agrees with the value found in 
[13] for area-preserving maps. 



N 


Ai -0 


A 2 


A3 - /u/A 2 


A 4 - n/ti 


7 


1.5 x 10" lb 


-0.7074647886 


1.2 x 10" 4 


5.2 x 10" 5 


9 


1.5 x 10" lb 


-0.7068366347 


7.8 x lO" 6 


3.4 x 10" y 


11 


1.5 x 10" 16 


-0.7067965675 


1.7 x 10" 7 


7.2 x 10" 8 


13 


1.5 x 10~ w 


-0.7067957946 


2.0 x 10" 8 


8.7 x 10- y 


15 


1.5 x lO" 16 


-0.7067956886 


4.5 x 10" 11 


1.6 x 10" 17 



Table 2. Scaling parameters \j 



In order to visualize the critical scaling Sh, , we computed some trajectories for the 
reduced flow in the variables q and z = O • p, 

(q, z) = (u + d 2 h(q, z)Q , - Q ■ Vi%, zj) , (2.6) 

for the Hamiltonian Hi 15) . The reduced phase space is considered to be T 2 x M, that is, 
q is identified with q + lixH? . Figure 1 shows 14 orbits of the return map for a rectangle 
in the plane IT = {(sv,z) : s,z G K}, where v = (—0.899,1). The coordinates in the 
horizontal and vertical direction are v ■ q and z, respectively. In these coordinates, the 
scaling fixed point is at (0, Zqo), where z^ = fl-pao = 0.0478 — The curve C through this 
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point represents the intersection of the critical cu-torus with IT. We note that v is not an 
eigenvector of A. The vector v has been chosen in such a way that C has zero curvature 
at the point (0, Zqo). In order words, we used the flow to linearize Sh„ approximately, in 
the angle variable. Figure 2, which will be discussed later, indicates that for most other 
choices of v, C would have a nonzero curvature at q = 0. As a result, the vertical scaling 
off the line v ■ q = would be dominated by A| instead of A3 . The points in Figure 1, 
that lie on the symmetry line v ■ q = and are surrounded by circles, belong to periodic 
orbits with periods $5 = 13/8, $6 = 21/13, .... If a similar picture holds for the first 
few #fc-orbits of H* , then the corresponding points (0, Zk) have to accumulate at (0, Zqq), 
with an asymptotic ratio (zk+i — z <x )/(zk — 2<x>) —> A3 = —0.326 . . .. And in the direction 
(v,0), the orbits have to scale with an asymptotic ratio equal to A2 . These properties 
follow from the fact that H* oSh* is a constant multiple of H* , which implies (since Sh„ is 
canonical modulo a scaling of the actions) that if an orbit for H* is periodic with frequency 
vector v, then its image under Sh* is another periodic orbit for H* , with frequency vector 
T~ x v. The observed accumulation of -orbits also suggests that the point (0,Poo) lies 
on an invariant u;-torus for H* , and that this torus is mapped onto itself by Sh* • For 
corresponding results on area-preserving map, we refer to [10-13]. 

The results obtained here and in [5] suggest that, as far as the breakup of smooth 
cu-tori is concerned, the unstable manifold of the transformation 1Z described after (1.3) 
represents the "typical" one-parameter family of Hamiltonians. On the other hand, these 
Hamiltonians are all of the form (2.1). Their cu-tori at different energies have to break up 
simultaneously, since they are just translates of each other. In order to see in what sense 
this may be typical, consider some "arbitrary" Hamiltonian H(q,p) that is even in q and of 
the form (1.1). Suppose that H allows us to define a differentiable one-parameter family of 
area-preserving maps E 1— > ^ e> by restricting the flow for H to surfaces H(q,p) = E and 
choosing an appropriate Poincare section. Assume now that H has a critical cu-torus for an 
energy E^ , and that the family E 1— > W e intersects the critical manifold W for the area- 
preserving renormalization group fixed point transversally, at E = E^ . Then, according to 
the renormalization group picture for area-preserving maps, there is an asymptotic scaling 
that relates the orbits of ^e to orbits of ^ f(E) 7 where f{E QO + e) = E^ + 62s + 0(s 2 ). 
Since these orbits all stem from a single Hamiltonian H , just considered on different energy 
surfaces, one would expect to find an asymptotic scaling rate S^ 1 = 0.614 ... in the phase 
space of H, in a direction transversal to the surface H(q,p) = E^ . But our results yield 
a scaling rate A4 = 0.142... in this direction. This suggests that the family E 1— > ^e 
cannot be transversal to W. In addition, if E 1— > ^e and W are smooth, then the fact that 
62 1 is not an integer power of A4, seems to imply that the distance between \& e and its 
projection onto W tends to zero faster than any power of E — E^ , as E approaches E^ . 
In the analytic case, this would lead to the conclusion that H has either no critical torus, 
or else an entire family of critical tori. We have not done any numerical investigations to 
confirm such behavior, or to find a counter-example. 

We conclude with some observations concerning the critical cu-torus for the fixed point 
H* , on the energy surface containing the point (0,Poo). A numerical approximation of this 
torus, in the coordinates of the reduced flow (2.6), is shown in Figure 2. These data suggest 
e.g. that the critical torus of H* is the graph of a differentiable function p = G(q), with 
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G(0) = Poo , and that is a "saddle point" of the corresponding function z = Q, • G(q). 
Notice that by symmetry, Q ■ G has to have a critical point at q = 0. The fact that the 
unstable eigenvector of DSn^iO^Pco) is tangent to the flow at (0,Poo), suggests that the 
orbit of the point (0,poo) coincides with the unstable manifold of Sh, at this point, and 
this in turn suggests that the critical w-torus is an attractor for the scaling map Sh* ■ In 
order to obtain information about the regularity of this torus, it is useful to look at periodic 
points of Sh„ ■ Consider e.g. the fixed point (0,p*). Our earlier discussion of Sh„ suggests 
that the degree of differentiability of the critical o;-torus at this point is no more than 
In | A3 1 / In I A2 1 = 3.229 . . .. The same holds of course for every point on the orbit of (0,p*) 
under the flow defined by H* . But there are other (dense sets of) points where the degree 
of differentiability is expected to be less than 2. Consider e.g. the set of points on T 2 
that are invariant under the reflection q 1— > —q. This set consists of 70 = (0, 0), and three 
other points 71,72,73, whose coordinates are both either or tt, modulo 2n. Notice that, 
according to Figure 2, these are precisely the critical points of O • G. It is easy to see that 
Uu, leaves the planes Tj = {7^} x IR 2 invariant. Thus, since T acts as a permutation on 
{71, 72, 73}, the corresponding planes Ti are permuted by Sh* ■ Assuming that the critical 
cu-torus is invariant under Sh* , this implies that {G(7i), £(72), ^(73)} is an orbit of period 
3 for Sh* • An estimate on the degree of differentiability of the critical torus at G(7i) can 
now be obtained by linearizing S^j at this point, and then computing the analogue of the 
ratio In I A3 1 /In I A2 1 . Based on observations for area-preserving maps [11,12], we expect 
the value of this ratio to be 1.790 . . ., but we did not verify this numerically. We note that 
for every periodic orbit of T on the torus T 2 , there is a corresponding periodic orbit of the 
transformation 3?. This is a consequence of the following property of Assuming that 
the equation (1.6) has a unique solution, it is easy to verify that SRo ,/ 7 = J T ~\ 1 o5R, where 
J 1 denotes the rotation H 1— > H(. — 7, .). 



3. Conclusion 

The renormalization group transformation 3? introduced in this paper can be used to study, 
and possibly explain, the breakup of golden invariant tori in Hamiltonian systems with two 
degrees of freedom. We find an approximate numerical fixed point for 3?, whose critical 
indices and scaling parameters are in excellent agreement with the corresponding values for 
area-preserving maps. The same methods are expected to work for other reduced quadratic 
irrationals, with similar results. Two distinguishing features of our transformation are 
that it applies to Hamiltonians of a general type, and that it seems to be analyticity 
improving. On a conceptual level, this indicates that all irrelevant degrees of freedom are 
being eliminated under the iteration of 3?. Numerically, the analyticity improving property 
allows for high accuracy, as truncation errors decrease exponentially in the degree of the 
truncation. Our results not only indicate that Hamiltonians with critical invariant w-tori 
are asymptotically self-similar, but they also provide detailed geometric information about 
the relevant scalings and normal forms. Further investigations will be necessary in order 
to determine the implications of these findings on the mechanisms underlying the breakup 
of invariant tori. Another project for future research is to carry out a computer-assisted 
proof for some of the results presented in this paper. 
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